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Presently only 30% of short gamma ray bursts (SGRBs) have accurate redshifts, and this sample 
is highly biased by the limited sensitivity of Swift to detect SGRBs. We account for the dom- 
inant biases to calculate a realistic SGRB rate density out to z = 0.5 using the Swift sample of 
peak fluxes, redshifts, and those SGRBs with a beaming angle constraint from X-ray/optical ob- 
servations. Assuming a significant fraction of binary neutron star mergers produce SGRBs, we 
calculate lower and upper detection rate limits of (1-180) per Yr by an advanced LIGO and Virgo 
coincidence search. Our detection rate is compatible with extrapolations using Galactic pulsar 
observations and population synthesis. 
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1. Introduction 

Binary neutron star mergers (NS-NS) are favoured as the progenitors for short gamma ray 
bursts (SGRBs), based on the association of some SGRBs with an older stellar population, and 
host galaxy types [|l|, f||. Kicks imparted to NSs at birth will produce velocities of several hundred 
km s _1 , implying that binary inspiraling systems may occur far from their site of origin. Fong, 
Berger & Fox [||] using Hubble Space Telescope observations to measure SGRB-galaxy offsets, 
find the offset distribution compares favourably with the predicted distribution for NS-NS binaries. 

Within 5 years, co-ordinated gamma-ray, X-ray, optical and gravitational-wave observations 
may allow the strong gravity regime of the central engine of compact object mergers to be probed. 
Such 'multi-messenger' observations provide the opportunity to probe these events across a vast en- 
ergy spectrum, and to constrain the progenitor populations of SGRBs. Furthermore, co-ordinated 
optical and gravitational-wave searches may play an important role in confirming the first direct 
gravitational-wave observations of compact object mergers [Q]. It is becoming increasingly impor- 
tant to constrain the rate of compact object mergers and their proposed optical counterparts in the 
context of up-coming gravitational-wave searches. 

We calculate a beaming-corrected SGRB rate density using the Swift sample of SGRB peak 
fluxes, redshifts and inferred beaming angles from X-ray observations. For our GRB selection 
criteria we use the Jochen Greiner catalogue of localized GRBs (see Table 1) and select bursts 
indicated as short that have reliable redshifts up to 2012 April. We avoid using a SGRB luminosity 
function, models for progenitor rate evolution, and a beaming angle distribution, all of which have 
large uncertainties. Instead, we focus on observed and measured parameters that take into account 
selection effects that modify Swiff s detection sensitivity to SGRBs. Finally, we use our SGRB rate 
density estimates to infer a detection rate of binary NS mergers by advanced LIGO (aLIGO) and 
Virgo (AdV) interferometers. Despite the poor statistics, this approach gives meaningful results 
and can be followed up when a larger sample of SGRB observations becomes available. 

2. Model 

The intrinsic SGRB rate is calculated taking into account the following observational bias 
effects using the sample in Table 1. GRB051221A is the only SGRB with a directly measured 
beaming angle in the sample (only a lower limit for the beaming angle was obtained for GRB 
070724A-see Coward et al. [Q]). We evaluated an upper SGRB rate density by using the smallest 
observed beaming angle in the sample (7 deg) and a lower rate by assuming isotropic emission. 

The low energy detection bandwidth of Swift (15-150 keV) in comparison with BATSE's 50- 
300 keV results in a bias against SGRBs which typically have harder emissions. Secondly, the Swift 
detection threshold is not simply defined by the detector sensitivity, but by a complex triggering 
algorithm Both these effects manifest as a bias against Swift detecting SGRBs; i.e. a smaller 
proportion of bursts has been detected by Swift (~10%) [gj. The latter effect results from the 
detection process employed by BAT, the Swift coded-aperture mask y-ray detector. In addition to 
requiring an increased photon count rate above background (the sole triggering criterion used for 
BATSE), BAT employs a second stage in which an image is formed by accumulating counts for up 
to 26 s [JJ. 
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Table 1: SGRB peak fluxes, Tgo and redshifts taken from the Swift online catalogue and 
http://www.mpe. mpg.de/^jcg/grbgen. html used to calculate Poisson rates. We use the 20-ms peak pho- 
ton fluxes from the BAT2 catalogue where possible - those marked by * are 1 -s peak photon fluxes. T We set 
a lower limit on the jet opening angle from the time when the Swift XRT monitoring stopped. "The proposed 
host galaxy at z — 0.1218 for GRB 080905A [|^| is a strong outlier to the Yonetoku relation [[llj] and a 
redshift z > 0.8 would make it consistent [ 12 1. See [E[ for expanded table with references. 



We attempt to crudely correct for these biases by using the observed SGRB rate from BATSE 
as a rate calibration for the Swift SGRB rate. Because BATSE operated at different energy thresh- 
olds and trigger sensitivities wn, for consistency we take all BATSE SGRBs with 64-ms peak flux 
when the trigger threshold was set to 5.5 a in the 50-300 keV energy range (total live operation 
time of 3.5 years). This yields 32 SGRBs sr _1 yr -1 , assuming an effective BATSE FoV of % sr [Qj. 

Because SGRBs occur over a short duration, it is more difficult (compared to long bursts) 
to produce a significant signal above background. Hence instead of using the theoretical BAT 
sensitivity of Fu m = 0.4 phs -1 cm~ 2 we employ a flux limit of 1.5 phs -1 cm~ 2 , using the 
smallest 20-ms peak flux from the SGRB data. 

Taking into account the sensitivity reduction and k-correction (see [Q] for a derivation of the 
k-correction), the SGRB all-sky rate can be inferred from the flux-limited SGRB sample in Table 
1. We calculate the maximum distance d max , with corresponding redshift z max , that a burst at 
luminosity distance di could be detected given Swift's sensitivity fLim (see for the derivation of 

Zmax)- 

The corresponding maximum SGRB detection volume for each burst is defined as 

V max = / ~rdz, (2.1) 
jo dz 

where the volume element, dV /dz, and luminosity distance, di(z), are calculated using a flat- A 
cosmology with Hq = 11 km s _1 Mpc -1 , Q.m = 0.3 and Q.^ = 0.7. 

To calculate the intrinsic rate of SGRBs requires accounting for the beaming angle, 0j, of the 



jetted burst. Equation ( |2.2| ) expresses the beaming factor used to correct for the unobserved SGRBs 



that are not detected because the jet is misaligned with the detector: 

B(0 j ) = [l-cos(e j )]- 1 . (2.2) 
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Table 2: The beaming-corrected SGRB rate densities with Poisson uncertainties using the observed con- 
straints on 0j, and scaled by the probability P,(7 go; p L ), except GRB 100816A, which uses just the T<)q. Lower 
rate estimates assume isotropic emission, and upper rates use the observed beaming angle constraints shown 
in Table 1, or the smallest observed beaming angle in the sample^, 0j » 7". " Because of the importance of 
GRB 080905 for the rate density, and its uncertainty (see Table 1 . for caveats), we also calculate total rates 
excluding this burst i.e. (3+^ - 500+J^) Gpc~ 3 yr -1 . 



To account for the fact that only a fraction of observed SGRB have measured redshifts, we 
scale the rate density by the ratio of Swift bursts with redshift to those without redshifts, F r m 8/39. 
The time span encompassing all observations, T ~ 6 yrs, is defined by the start of Swift observations 
to the time of the most recent SGRB in the sample and we account for the fractional sky coverage 
of Swift, £2 « 0.17. To account for Swift's reduced sensitivity for detecting SGRBs relative to 
BATSE, we use the ratio of the BATSE to Swift SGRB detection rate, which we approximate as 
Rb/s = 6.7. This converts the Swift SGRB rate to an intrinsic SGRB rate. We point out that this 
correction applies only to SGRB (see [Q] for an analysis including SGRBs with extended emission). 



Finally, we use the probabilities, P,(r 9o; p L ), that each burst is a non-collapsar [|13j, |14| ], and scale the 
rate density by these probabilities. We note that these probabilities do not affect the rate densities 
significantly. 

Combining all detection parameters yields the Poisson SGRB rate density of a single (/th) 
burst, and the total rate density for n bursts: 

2 i 111 

#SGRB = Yjw. FToV'^'^ft)' < - 2-3 ' ) 

i v i(max) r r 1 " 

3. Results and summary 

We use the above SGRB rate density to infer a detection rate of binary NS mergers by advanced 
gravitational-wave interferometers. For aLIGO and AdV interferometer sensitivities, the horizon 
distance, Dp, (all sky locations and orientation averaged over) for optimal detection of a NS-NS 
merger in a coincidence search is about 340 Mpc (cosmological redshift not included). For a direct 



comparison with [15], the detection rate is computed for a single interferometer, with D/, = 197 



Mpc, and for optimal detection in a three detector coincidence search (see figure nh. Given the 
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< 10 10 20 30 40 50 60 70 80 90 
SGRB beaming angle (degrees) 

Figure 1: The lower curve plots equation 3.1, the detection rate of binary NS mergers by a single aLIGO 
interferometer (D/, = 197 Mpc) as a function of SGRB beaming angle, using /?low = 8 Gpc~ 3 yr _1 (see 
Table 2). The upper curve assumes a three detector coincidence search with D/, — 340 Mpc. Both horizon 
distances used to calculate detection rates are angle averaged over all binary orientations. The dashed curve 
is the same as the upper curve but using the upper rate excluding GRB 080905 (see Table 1 . for caveats). 



uncertainty in the beaming angle distribution, we define a detection rate as a function of 6>j using 
the SGRB lower rate estimate, i.e. from 0j = 90°, scaled by B(dj) and the Euclidean volume: 

R(0- 3 ) = —DlR^BiGj). (3.1) 

For a realistic SGRB beaming angle range (0j = 7 — 30) degrees, the corresponding detection rate 
by a coincidence search can range from (10 — 180) detections per year. 

The binary NS gravitational-wave detection rate estimates are based on calculating an intrinsic 
SGRB rate density using Swift localized bursts, taking into account dominant selection effects. This 
approach, based on observational data is very different from that based on Galactic binary pulsar 



observations and modelled population synthesis. In the latter, [ |15[ ] use the observed Galactic binary 
pulsar population and extrapolate a NS merger rate density out to the aLIGO and AdV detection 
horizon. Conversely, our approach avoids this extrapolation because it is essentially an observed 
rate extending well beyond the average sensitivity distances of the upcoming gravitational-wave 
searches for compact binaries (about 300 Mpc or z = 0.07 for aLIGO and AdV interferometers). 

In conclusion, the upcoming gravitational-wave detection era will be fundamental for resolv- 
ing the SGRB-binary NS merger connection, since an unequivocal association between SGRBs and 
binary NS mergers will only be possible via coincident gravitational-wave and electromagnetic ob- 
servations. Ultimately, a comparison beteween the SGRB rate density and the gravitational-wave 
detection rate will help constrain the fraction of binary NS mergers that give rise to SGRBs and the 
SGRB beaming angle distribution. 
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